Dynamics of ripple formation in sputter erosion: nonlinear phenomena 
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Many morphological features of sputter eroded surfaces are determined by the balance between 
ion induced linear instability and surface diffusion. However, the impact of the nonlinear terms on 
the morphology is less understood. We demonstrate that while at short times ripple formation is 
described by the linear theory, after a characteristic time the nonlinear terms determine the surface 
morphology by either destroying the ripples, or generating a new rotated ripple structure. We show 
that the morphological transitions induced by the nonlinear effects can be detected by monitoring 
the surface width and the erosion velocity. 
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The morphological evolution of ion sputtered surfaces 
has generated much experimental and theoretical inter- 
est in recent years. As a result, there is extensive evi- 
dence that ion bombardment can result in ordered sur- 
face ripples, or lead to kinetic roughening, depending on 
the experimental conditions. These experimental results, 
which cover amorphous and crystalline materials (SiC>2 
|1|), and both metals (Ag and semiconductors (Ge 
[3 1, Si j§||]), have motivated extensive theoretical work 
aiming to uncover the mechanism responsible for ripple 
formation and kinetic roughening. A particularly suc- 
cessful model has been proposed by Bradley and Harper 
(BH) ||, in which the height h(x, y, t) of the eroded sur- 
face is described by the linear equation 



+ ^(d x h) 2 + ^-(d v h) 2 + V (x,y,t), 



(3) 



d t h 



v x d x h + Vyd y h 



Kd 4 h, 



(1) 



where K is the surface diffusion constant and the coeffi- 
cients v x and v y are induced by the erosion process such 
that v y < and v x can change sign as the angle of the in- 
cidence of the ions is varied. The balance of the unstable 
erosion term (— \v\d 2 h) and the smoothening surface dif- 
fusion term (—Kd 4 h) generates ripples with wavelength 



U = 2tt^2KJ\u~\ 



(2) 



where i refers to the direction (x or y) along which the 
associated Vi (y x or v y ) is the largest. While successful in 
predicting the ripple wavelength and orientation , this 
linear theory cannot explain a number of experimental 
features, such as the saturation of the ripple amplitude 
^|-|l0|, the observation of rotated ripples Jl|], and the 
appearance of kinetic roughening Jl^,[lj| ■ Recently it has 
been proposed ||l4|| that the inclusion of nonlinear terms 
and noise (both of which were derived from Sigmund's 
theory of sputtering p5| ) can cure these shortcomings. 
Consequently, Eq.(l) has to be replaced by the noisy 
Kuramoto-Sivashinsky equation (KS) equation |l6| 



where rj{x, y, t) is an uncorrelated white noise with zero 
mean, mimicking the randomness resulting from the 
stochastic nature of ion arrival to the surface. All coef- 
ficients in Eq.(3) have been determined in terms of the 
experimental parameters, such as the ion flux, angle of 
incidence, ion penetration depth, and substrate temper- 
ature 0,0. 

While it is expected that the nonlinear theory incor- 
porates most features of ripple formation and kinetic 
roughening, the morphological and dynamical features 
of the surface described by it are known only in certain 
special cases. For example, when the nonlinear terms 
and the noise are neglected (X x — \ y — 0, rj = 0), 
Eq.(3) reduces to the linear theory (1), and predicts rip- 
ple formation. It is also known that the isotropic KS 



equation (y x 



Vy < 0, K x 



K y = K xy /2, and A^ = X y ) 
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asymptotically (for large time and length scales) predicts 
kinetic roughening, with exponents similar to that seen 
experimentally in ion sputtering [ jl2| . For positive v x 
and v y , Eq.(3) reduces to the anisotropic KPZ equation 
p8| , whose scaling behavior is controlled by the sign of 
Ax • A y jl9| . Finally, recent integration by Rost and Krug 
p0| of the noiseless version of Eq.(3) provided evidence 
that when \ x ■ X y < 0, new ripples, unaccounted for by 
the linear theory, appear and their direction is rotated 
with respect to the ion direction j2(J. However, it is not 
known if this rotated phase survives in the presence of 
noise, rj. These special cases illustrate the complexity of 
the morphological evolution predicted by Eq.(3). To be 
able to make specific predictions on the morphology of 
ion-sputteted surfaces, we need to gain a full understand- 
ing of the behavior predicted by (3), going beyond the 
special cases, that are often experimentally irrelevant. 
In this paper, we integrate numerically Eq.(3), aiming to 
uncover the dynamics and the morphology of the surfaces 



for different values of the parameters. We demonstrate 
a clear separation of the linear and nonlinear behavior. 
For short erosion times, the ripple development and ori- 
entation follow the predictions of the linear theory of 
BH. However, after a well defined crossover time, which 
depends on the coefficients of Eq.(3), the surface mor- 
phology is determined by the nonlinear terms. We find 
that when X x ■ X y > the nonlinear terms destroy the 
ripple morphology. However, when X x ■ X y < 0, they 
result in a long and apparently rough transient regime, 
followed by a new morphology of rotated ripples, as seen 
in the noiseless KS equation. We show that these mor- 
phological transitions can be detected by monitoring the 
surface width or erosion velocity, quantities that can be 
measured more easily in situ. Finally, we discuss the 
impact of our result on current and future experimental 
work. 

The direct numerical integration is carried out by dis- 
cretizing the continuum equation of (3), using the stan- 
dard discretization techniques . We choose a tempo- 
ral increment At = 0.01 and impose periodic boundary 
conditions h(x, y, t) = h(x + L, y, t) = h(x,y + L, t) where 
L x L is the size of the substrate. We choose the noise 
to be uniformly distributed between [-1/2,1/2], and per- 
turb the initial flat configuration with the noise. Since 
the sign of the nonlinear terms plays a significant role in 
defining the surface morphology, we discuss separately 
the X x ■ X y > and X x ■ X y < cases. 

X x ■ X y > — A general feature of sytems such as 
Eq.(3) is that the nonlinear terms do not affect the sur- 
face morphology or dynamics until a crossover time r has 
been reached. Thus, we expect that for early times, i.e. 
for t < t, the surface morphology and dynamics is prop- 
erly described by the linear theory To demonstrate this 
separation of the linear and nonlinear regimes, in Fig. 1 
we show the time dependence of the surface width defined 
as W 2 (L, t) = J2 X y h 2 (x, y, t) — h 2 and of the mean 
height h—j^ ^ x h(x, y, t). We find that for t < t, the 
width W increases exponentially while the mean height 
stays constant at h = 0. Indeed, both of these findings 
are consistent with the predictions of the linear theory: 
W, being proportional to the ripple amplitude, according 
to Eq.(l) increases as W ~ exp(i/t/£ 2 ), and the linear 
terms do not change the mean height of the surface. 
Furthermore, inspecting the surface morphology, we find 
that in this regime the ripple wavelength and orientation 
are also correctly described by the linear theory. For ex- 
ample, for the parameters v x = —0.0001, v y = —0.6169, 
K x = K y = K xy /2 = 2, and X x = X y = -0.001, ac- 
cording to (2), the ripple wavelength along the y axis 
is i y w 16, and along the cc-axis is £ x w 1257. Since 
the dominant wavelength is determined by the growth 
rate ~ ex-p(vt / £ 2 ) , the smaller wavelength, i.e. £ x , will 
dominate. As Fig. 2a shows, for a system of size 64 x 64 
we observe four ripples aligned along the x-axis, in agree- 



ment with the previous prediction. As a second example 
we consider the case, v x = —1.2337, v y = —0.0001, 
K x = K y = K xy /2 = 1, and X x = X y = -0.001, for 
which we expect ripples of wavelength £ x » 8, smaller 
than £ y w 889. As Fig. 2b shows, in this case we observe 
eight ripples aligned along the y-axis. 

While the early time behavior is correctly predicted by 
the linear theory, beyond the crossover time r the non- 
linear terms become effective. One of the most striking 
consequence of these terms is that the surface width sta- 
bilizes rather abruptly (see Fig. 1). Furthermore, the rip- 
ple pattern generated in the linear regime disappears, and 
the surface exhibits kinetic roughening. A typical sur- 
face morphology, demonstrating the absence of ripples, 
is shown in Fig. 2c. The crossover time r from the linear 
to the nonlinear behavior can be estimated by comparing 
the strength of the linear term with that of the nonlin- 
ear term. Let the typical height at the crossover time r 
be Wq = 1/W 2 (L,t). Then, from the linear equation we 
obtain Wq ~ exp(z/r/£ 2 ), while from d t h ~ X(dh) 2 we es- 
timate Wq/t ~ XW 2 /£ 2 . Combining these two relations 
we obtain 

t ~ {K/v 2 )\n{v/X). (4) 

In this expression, v, K and A, refer to the direction 
perpendicular to the ripple orientation. The predicted 
A-dependence of r is confirmed in the inset of Fig. la. 
An another quantity that refects the transition from the 
linear to the nonlinear regime is the erosion velocity 
v = dth. The main contribution to the erosion veloc- 
ity comes from a constant erosion rate — vq, that has 
been omitted from (1) and (3), since it does not affect 
the surface morphology ]qjl4j]- However, in addition to 
vq, the mean height is also modified by the nonlinear 
terms, X x (d x h) 2 and X y (d y h) 2 . In the following for sim- 
plicity, we neglect the vo term, since its value does not 
depend on the surface morphology, and it is constant 
throughout the erosion process. The nonlinear terms 
act to decrease the mean height in the case of X x < 
and X y < 0. We can estimate the surface velocity as 
v ~ XWg/£ 2 - v z /{KX) using W - vjX. This depen- 
dence of v on A is consistent with the numerical results, 
shown in the inset of Fig. lb. In this regime (t > r) the 
surface exhibits kinetic roughening, i.e. the surface width 
should increase either logarithmically (when A = 0) or as 
a power law (when A ^ 0) fl9|p2[ . However, compared 
with the exponential increase in the early regime (t < r) , 
this dependence is hardly observable. The simulation 
times required to investigate the asymptotic scalings of 
W with t are currently prohibitive. 

X x ■ X y < — As Fig. 3a shows, we again observe a 
separation of the linear and nonlinear regimes, however, 
we find that the morphology and the dynamics of the 
surface in the nonlinear regime is quite different from the 
case X x ■ X y > 0. In regime I, for early times (t < r), 
the surface forms ripples (see Fig. 4a), whose wavelength 
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and orientation is correctly described by the linear the- 
ory. After the first crossover time r, given by Eq. (4), 
the surface width is stabilized, and the ripples disappear, 
as shown in Fig. 4b. After r, the system enters a rather 
long transient regime, that we call the regime II. Here, 
the surface is rough, and no apparent spatial order is 
present. We often observe the development of individual 
ripples, but they soon disappear, and no long-range order 
is present in the system. However, at a second crossover 
time T2, a new ripple structure suddenly forms, as shown 
in Fig. 4c, in which the ripples are stable and rotated with 
an angle 9 C to the x direction. The angle 9 C has the value 
6 C = tan -1 yj—\ x /\ y , (or tan -1 y/-X y /X x ) @, which 
can be calculated by moving to a rotated frame of coordi- 
nates that vanishes the nonlinear term in the transverse 
direction. Indeed, as Fig. 4c shows, the observed angle is 
in excellent agreement with 6 C = tan (1/2) for X x = 1 
and A y = —4. We also find that the time the system 
spends in regime II fluctuates from system to system, 
thus T2 has a wide distribution. The transitions between 
the three regimes can be detected by monitoring the sur- 
face width W (Fig. 3a): in regime I, the width increases 
exponentially, as predicted by the linear theory; it is ap- 
proximately constant but hightly fluctuating in regime 
II, and suddenly increases and stabilizes in regime III. 
Note that the amplitude of the rotated ripples in regime 
III is much larger than in regime II, a rather attractive 
feature for possible applications as patterned templates 
for various microelectronic applications. 

The demonstrated morphological transitions generate 
an anomalous behavior in h as well. As Fig. 3a shows, 
the mean height is zero in the linear regime, increases as 
the ripples are destroyed in regime II, and decreases with 
a constant velocity in regime III. In order to understand 
this complex behavior, we consider a specific example, for 
which the surface morphologies are shown in Fig. 4. For 
this parameter set, ripples are aligned along the y-axis in 
the region I, because l x <C l y . Thus, the contribution of 
(d x h) 2 is much larger than that of (d y h) 2 , even though 
|A X | < |Aj,|, and the surface height increases due to the 
term X x (d x h) 2 with A x > in regime II. However, as the 
ripples are destroyed by the nonlinear effects, the contri- 
bution of {dyh) 2 term increases, and eventually X y (d y h) 2 
becomes larger than X x (d x h) 2 , forcing the mean height 
to decrease because X y < 0. The velocity in regime III is 
determined by the nonlinear coefficient in the direction 
along the ripples, which reduces to Aa, + X y after the co- 
ordinate transformation to the rotated ripple direction. 
This prediction is in good agreement with the results of 
Fig. 4, which demonstrates that v ~ 1/(A X + X y ). 

The clear separation of the linear and the nonlinear 
behavior, that holds for both signs of X x ■ X y , has a direct 
impact on the experimental observations. Numerous 
experiments have observed the development of ripples 
whose wavelength and orientation is in good agreement 
with the prediction of the linear theory |l]-|!,i|,En] ■ Based 



on our results, we expect that these experiments were in 
the t < t regime, where indeed the linear theory fully 
describes the system. However, recent results have pro- 
vided detailed experimental evidence of ripple amplitude 
stabilization [p|-^0[ , a clear sign of the presence of nonlin- 
ear effects. Furthermore, it was found that the different 
W versus t curves can be collapsed by rescaling time 
with a factor v 2 jK and amplitude with \Jv/2K J9|. In- 
deed, this is in excellent agreement with our prediction, 
Eq. (4) . Finally, since the values of v and A can be tuned 
by changing the ion energy and the angle of incidence, 
and K can be tuned with the temperature, the values 
of r and T2 can be changed continuously, and thus our 
predictions on the morphological transition between the 
linear and nonlinear regimes could be tested experimen- 
tally. Furthermore, the detailed morphological evolution 
uncovered here, combined with earlier calculations that 
connect the coefficients in Eq.(3) to the numerical val- 
ues of the parameters describing the ion-bombardment 
process j(||l4|,[l7| , offer a detailed roadway that can guide 
further experiments and facilitate the use of ion sputter- 
ing for surface patternings. 
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FIG. 1. Time evolution of (a) the surface width W 2 and 
(b) the mean height h for the parameters v x = —0.0001, 
v y — —0.6169, K x = K y = K xy /2 — 2. The different curves 
correspond to different values of X x — X y — A. In (a), from top 



to bottom, the curves correspond to A = 0, —10 



-10" 



— 10~ 3 , — 10~ 2 , and — 10 _1 , respectively. In (b), from bot- 
tom to top, they correspond to A = — 10~ J , — 10~ 4 , — 10~ 3 , 
— 10 -2 , and — 10 _1 , respectively. Inset (a): The crossover 
time r, estimated from (a) is shown as a function of ln|A|. 
Inset (b): Plot of ln \v\ versus ln|A|. The dotted line has a 
slope ps —1.07, implying v ~ 1/A. 




(a) (b) (c) 

FIG. 2. Grey scale plot of a surface of size 64 x 64 with 
the parameter sets: (a) v x = —0.0001, v v — —0.6169, 



K 



K v 



-0.001 at t = 2x 10 4 



xy /2 = l and 



and X x 

(b) v x = -1.2337, vy = -0.0001, K x = K y = K, 
X x = X y = -0.001 at t = 2 x 10 3 . Eq. (2) predicts the wave- 
length £ y ?s 16 for (a), and 4 ~ 8 for (b). (c) The surface 
configuration at time ( = 5x 10 5 for the parameter set used 
in (a). 



2000 


-4000 
-8000 
-12000 



- (I) 


(II) 


(ni) 






- — h 



















- 10 4 



10 3 



10 5 
Time 



io 2 



- 10° 



JU- 



KI" 

io- 1 
\v\ 



10 7 



IO" 1 10" 
\Xx+Xy\ 



(a) 0» 
FIG. 3. (a) Time evolution of the mean height h (dashed, 

left linear scale) and the surface width (solid, right logarith- 
mic scale) for the parameters, v x = —0.6169, v y = —0.01, 
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K xy /2 — 2, A a 



1 and A,, 



-4. The dotted 



lines seperate the three regimes discussed in the text, (b) 
The dependence of |u| on the nonlinear terms \X X + X y \ for 
the same parameters used in (a). The dotted line has a slope 
ps —1.02, implying v ~ l/(Ai + A2). 
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(a) (b) (c) 

FIG. 4. Grey scale plot of a surface of size 64 x 64 with the 
parameter set used in Fig. 3, showing the surface morpholo- 
gies at t — 10 4 in regime I (a), t = 2 x 10 5 in regime II (b), 
and t = 10 7 in regime III (c). 
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